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Abstract 

We study the problem of relaxing a large cosmological constant in the astrophysical domain 
through a dynamical mechanism based on a modified action of gravity previously considered by 
us at the cosmological level. We solve the model in the Schwarzschild-de Sitter metric for large 
and small astrophysical scales, and address its physical interpretation by separately studying the 
Jordan's frame and Einstein's frame formulations of it. In particular, we determine the extremely 
weak strength of fifth forces in our model and show that they are virtually unobservable. Finally, we 
estimate the infiuence that the relaxation mechanism may have on pulling apart the values of the two 
gravitational potentials ^(r) and <l?(r) of the metric, as this implies a departure of the model from 
General Relativity and could eventually provide an observational test of the new framework at large 
astrophysical scales, e.g. through gravitational lensing. 

1 Introduction 

The old cosmological constant (CC) problem [1] is considered as one of the biggest puzzles in theoretical 
physics. Remarkably, it is only partially related to the question of what drives the current accelerated 
expansion |2ji3i|4j. In general the latter question is "answered" by the so-called "dark energy" (DE), which 
summarizes a large amount of different models and concepts |5J [51171 [HI IHJ [ID] • The simplest candidate for 
the DE is a tiny positive cosmological constant A, others are quintessence scalar fields, modified gravity, 
the existence of large scale ambiguities or even a misinterpretation of observations, just to mention a few 
of them. Whatever the true answer might be, the old CC problem still remains in its full glory because 
we cannot just put the huge vacuum energy predicted by every quantum field theory (QFT) under the 
rug, and we cannot just resolve the problem by mere fine tuning methods - cf. e.g. the recent references 
HB [H US] and [HQ In essence, the CC problem stems from the profound incompatibility between the 
observed DE density [3J [1] with the theoretical expectations [U [5J [3 [HI 13 HO] ■ One has to admit that 
theoretical physics is currently not able to predict the value of the vacuum energy density pA — A/ (SttG). 
In actual fact, this is not the main preoccupation, since after all QFT is unable to predict, say, the value 
of the electron mass. The real problem instead is quite another one, to wit: while renormalizable QFT 
does provide a fully consistent framework to accommodate the value of the electron mass (namely, one 
which is free from UV and IR ambiguities) , it is nevertheless unable to do the same (as far as we know) 
for the measured value of pA- As a result, we are left with various rough order of magnitude estimates 
of finite vacuum energy contributions, e.g. from phase transitions and inflation, in addition to naive and 
mostly unreliable calculations of quantum zero-point energies - see however |16l [171 [IH] and references 
therein. In this untidy state of affairs, there is at least one single thing of which we can say that we are 
fully convinced to be true: that the huge contributions to the CC come from theories that work at typical 

^See specially the detailed account of the fine tuning CC problem in Appendix B of |14 |. 



energies much larger than those observed in the present universe, and therefore any practical mechanism 
is bound to start working effectively and efficiently already at these large scales, and then persist in the 
job for a sufficiently large period of time until our days. 

A promising dynamical route along these lines, i.e. one which is free from fine tuning problems, is 
taken by models which intend to make a low-energy universe as ours feasible despite having a large CC in 
the total energy content. In Refs. |11|[T^[T51[H] we proposed some intriguing models of this kind, where 
the machinery at work was called the CC relaxation mechanism. The resulting cosmos was referred to 
as the relaxed universe. It is not surprising that these more innovative approaches require more invasive 
changes in the theory than just providing a source of late-time acceleration, as it occurs in the usual 
modified gravity models [19l|20lllll[22j|23l[211l25l|26]. While some of the latter are well motivated by 
observations, they usually presume the existence of an implicit (and extremely fine tuned) counterterm 
in the effective action just devised to cancel the big CC, such that the other sources (matter and the new 
terms of the modified gravitational action) can drive the acceleration on top of the essentially zeroed (by 
hand) vacuum energy. Hence, even the most viable modified gravity models in the literature are in urgent 
need of some kind of CC relaxation that can avoid the unacceptable fine tuning of the vacuum energy 
which is implicit in all of them. The enormous difficulty in providing some form of solution capable 
to alleviate this acute problem plaguing all the traditional (i.e. "late-time") models of modified gravity 
gives a strong motivation to invoke new and unconventional concepts, with the understanding that in the 
beginning they might have the status of toy models or prototype ideas just hinting at the final solution. 
A serious hint along these lines might just be to replace the "late-time" modification of gravity by an "all 
time" (i.e. perennial or everlasting) modification of it. 

In Ref. |14| we discussed a specific perennial modified form of gravity which provides an implementa- 
tion of the CC relaxation mechanism at the cosmological level using the FLRW metric. We found that 
the evolution of the universe is well behaved despite the presence of a huge initial vacuum energy density 
of order p\ ^ Mjr < Mp {Mp being the Planck mass), in fact much larger in magnitude than the energy 
density of matter and radiation, thanks to the dynamical compensation of p\ triggered by the expansion 
itself 1^ The fact that this relaxation mechanism must be active during most periods of cosmological evo- 
lution (actually, at all times after the infiationary period) shows the intriguing nature of this approach. 
But it also explains the difficulty of generalizing the mechanism to more complicated systems, specially 
those that unavoidably involve the interplay of gravity and the matter sources. In order to understand 
better the properties of the CC relaxation mechanism developed in [TJ , we investigate here its behavior 
on sub-horizon scales, actually on length scales fully in the astrophysical domain, e.g. in the Solar system 
environment. Apart from searching for adequate metric solutions in spherically symmetric backgrounds, 
and dwelling on the physical interpretation of our model in different gravity frames, we will estimate 
the strength of the possible fifth forces. Indeed, in many modified gravity theories these forces appear 
as a manifestation of extra gravitational degrees of freedom (d.o.f.). The purpose of this Letter is to 
obtain a qualitative understanding of all these effects. Hence, we restrict the discussion to the special 
case of static backgrounds. At the cosmological level this corresponds to the de Sitter solution we found 
in [TJ, which defines the border line between the quintessence-like and phantom-like behavior exhibited 
by our model in the future regime. Consistently, we expect that the static solutions we deal with do not 
represent completely stable configurations. However, since we have found in [14J that the cosmological 
future can be very close to the de Sitter cosmos at the present time, we may nevertheless obtain useful 
information from the study of the static case without facing the full time-dependent situation. Finally, 
we discuss some numerical solutions leading to possible deviations with respect to General Relativity. 



2 The relaxation mechanism at astrophysical scales 



The action describing our scenario is given by |14| 
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with the standard Einstein-Hilbert term {l/2)MpR — i?/(167r Gjv), the matter action S'mat and the 
large vacuum energy term p^. The modifications to gravity are described by terms proportional to j3 



detailed confrontation of these kind of models with the modern cosmological data has been recently presented in 
Ref. 1151 . It is remarkable that they are able to fit the expansion history of the universe in a comparable way to the 
concordance ACDM model. 
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and (3q, and the function B = B{R,G) is given in terms of the Ricci scalar R and the Gaufi-Bonnet 
invariant G = - ARabR'''' + RabcdR"'""^ ■ 

The corresponding gravitational field equations follow from the functional variation of the action (jlj. 
After a straightforward calculation one finds: 

Rab - ^9abR = SttGn [Tab + Pa 9ab + 2i?afc] , (2) 

where Tab is the ordinary energy-momentum tensor, and Eab is the extra tensor that appears as a 
consequence of the modified gravity terms. Explicitly, 

Eab = -gab{j^ + ^j+{Rab + VaVb-gabO)l-—R+ — 

p \ / pn r. 



+ □^-^j_V™V„(^^j|. (3) 

We see that this tensor follows solely from the gravity modifications in Q. If we set (3 = f3Q = in 
it, then Eab = and Q boils down to the standard Einstein's equations with non- vanishing vacuum 
energy density p\, i.e. Rab — {^/2)gabR — —SttG (Tab + Pa 9ab)- In the general case, the field equations 
are rather cumbersome and we will have to consider different simplified situations. 

A working ansatz for B{R, G) is given in Ref.[T3], where it is then specialized to the FLRW background 
and hence valid for cosmological considerations. In this domain, and for all epochs after the radiation 
epoch, it takes the simplified form 

B{R,G) = ^R^ + ^G, (4) 

which is the precise combination that reproduces the correct radiation and matter epochs through the 
cosmological evolution [T3]. Therefore, this expression also applies for astrophysical applications, although 
in this case it must be evaluated for a metric amenable to the new context, typically a spherically 
symmetric one (see below). 

On theoretical grounds, the initial CC density in the early universe should read p\ ^ Mj^ (where Mx 
can be some GUT scale near Mp), and is to be taken much larger in magnitude than any other energy 
density. Therefore, the MpR and matter terms pmat in the action ([l]) are much smaller than p\. Even 
if we take p\ as the vacuum energy density of the Standard Model (SM) of Particle Physics, we have 
p\ ~ —(2 X 10^ GeV)'*, which induces a very large cosmological constant, some 56 orders of magnitude 
away from the measured value ~ 10~*^ GeV^. That (negative) contribution from the electro- weak 
phase transition is probably the most reliable number we can assume at the moment for the theoretical 
value of the vacuum energy, inasmuch as it represents the prediction of the SM as the most successful 
QFT to date. According to our earlier work [T5J [T31 [T3], the gravity modifications in ^ can be 
used to relax dynamically the large CC and leaving a low-curvature cosmological solution which is not 
dominated by p\. This is possible because the effective vacuum energy density in our F{R, G)-cosmology 
- and hence the quantity playing the role of DE in our framework - is not just the parameter p\^, but 
the full expression pAcff = Pa + 2i?o, as one can see from The effective quantity pAcff is dynamically 
enforced to be very small by the relaxation mechanism [14J. If we focus on the relevant terms of the action 
which are acted upon by this mechanism, the corresponding field equations read 

2Eab + p\9ab^0. (5) 

Solving this equation means that there is a dynamical cancelation of p\ by the /3,/30-terms in The 
remaining terms play no role for the relaxation because 



B R 



|pX I > Pmat- (6) 



The details of the relaxation mechanism and various numerical examples in the cosmological context have 
been discussed at length in the comprehensive paper [T3]. 

For typical values of the /3 parameters, the Pq/R term has little infiuence on the large scale cosmo- 
logical expansion (as we shall discuss below). Therefore, the standard problems of 1/R gravity (as e.g. 
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the incorrect description of the matter and radiation dominated epochs in these models |27l B5] ) do not 
apply in our case because at very large (cosmological) scales the leading term in the action ([l]) is not 1/i? 
but l/B, and we know that the latter works perfectly well since it is capable of relaxing the large CC 
and it correctly reproduces the standard radiation and matter dominated epochs [14]. 

On small length scales, instead, it is the 1/i? term of ^ that dominates, and therefore it will be this 
term that is responsible for the CC relaxation in the astrophysical domain. To better understand this 
behavior, let us study static spherically symmetric solutions described by the metric 

ds^ = a(r)df -TT^-r^ (dd^ + sin^ 9 dcf^) , (7) 
b(r) 

with two functions a(r), h{r) depending only on the radial coordinate r. The time and angle coordinates 
play no role in the following. As mentioned already in Sec. [T] static solutions on small scales correspond 
to the asymptotic de Sitter metric at the cosmological level separating the quintessence from the phantom 
future behavior. Consequently, the static setup considered here is not expected to be absolutely stable. 
However, it may be a sufficiently useful approximation for obtaining some insights into situations not far 
to the de Sitter future. According to the recent observations, which prefer a CC-like equation of state for 
the DE, lod ~ —1, this a reasonable assumption. 

In the general case, finding solutions of the relaxation equation (|5| for the astrophysical problems can 
get quite complicated, so further approximations are necessary. The simplest situation that we wish to 
consider is the de Sitter, spherically symmetric, space-time. It corresponds to the form given in Eq. ([t]) 
with 

a(r) = l-^, &(r) = a(r), (8) 

thus representing a 3-dimensional sphere of radius — y^S/A, which acts as a cosmological horizon. 
Such a space-time is an exact solution of (|5| for all r, as we can check. Let us indeed compute the extra 
tensor Eab for the de Sitter metric. First of all we find R = 12r~^ and B — lQ8r~'^, where we used Q 
to compute the latter. Finally, from ([3| we arrive at: 



- 1^108 + ^« 16 

The fact that Eat is constant confirms that the de Sitter space-time is an exact solution. Equation ([s]) 
can then be satisfied and so the expression (|9|) cancels against ^p\gab- As both terms with (3 and /3q 
are in general needed to cancel p\, they together need to produce the contribution of the sign opposite 
to p\. The straightforward way to achieve this for > is to have both /3 < and I3q < 0, whereas for 
p\ < to have both /3 > and /3q > 0. Finally, after imposing the boundary condition on the current 
value of A, this fixes the order of magnitude of these parameters]^ Of course the cancelation is not exact, 
leaving a very small remainder (the current value of p^), but the mechanism does not depend on it. Let 
us emphasize that this cancelation involves no fine tuning because rg is not a constant to be fixed by 
us, but a dynamical variable controlled by the relaxation mechanism. Whatever it be the starting value 
for p\, the relaxation mechanism chooses such that (jsj) is fulfilled. Recall that = ^^S/A = H^^ 
(with A = SttGn PAeff) is driven to the current value of order Hq^. The reason for this dynamical choice 
stems from the late time epoch of the universe evolution, in which the relaxation condition (viz. i? — ^ 0, 
without ever being exactly zero) enforces a very small value of H |14| . Since becomes then very large, 
we see from Eq.(pl that the expression Eqq becomes also very large, in fact as large as to essentially cancel 
against p\ ^ Ai/j- Notice that the first term of ^ behaves as ^ and hence dominates at very large 
distances; this is the term that emerges from the 1/S-invariant and which we used in the cosmological 
domain to insure a very small value for the measured CC |14| . As we can see from in this domain 
is dynamically driven to satisfy /3 rg/54 -I- p\ = p° <C p\, i.e. effectively Eq. ([s]). Therefore, 



Pa 



^ M 



(10) 



■'it should be clear that the values of /3 and need not be fine-tuned at all, we only have to fix their order of magnitude 
and sign - see equations ( |10| and l |16| below. The relaxation mechanism then selects re automatically such that |5]l is 
fulfilled, irrespective of the input value of pj^ 1141 
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where we have set \P\ = M^. Choosing M = ^/MxHq we get - H^'^ . For Mx ^ 10^'^"^'^GeV 
(GUT/string scale), we see that the required order of magnitude for Ai is in the ^meV range, charac- 
teristic of a Hght neutrino mass, and is also the energy scale of the CC: ttia = (pa.)^'''' ^ 10^'^ eV. 

The relaxation mechanism also works fairly well in the astrophysical domain. It only changes qualita- 
tively, as we shall discuss below, where now it is the 1/R term of the action ([ij that takes over. But the 
relaxation mechanism once more works appropriately to protect the astrophysical scales from unwanted 
vacuum effects. To start with, we note that the astrophysical domain cannot be described by just the 
spherically symmetric de Sitter metric ([8|; we need the Schwarzschild component too. Let us therefore 
consider the Schwarzschild-de Sitter (SdS) ansatz with non-zero Schwarzschild radius <C r^: 

= = (11) 

r 

In contrast to standard General Relativity, this metric is not a solution of our system because the 
Gaufi-Bonnet term in j3/B induces an r-dependence in the extra tensor Eab in ([s]), which cannot be 
compensated by the constant value of p\. However, for sufhciently small or large values of the radius, 
Eab is still approximately constant. To better understand this issue, let us look at the extra tensor 
following from this metric. In the astrophysical setup, the invariant functions in the action ([T|) take on 
the form 

R = 12r-^ , B = lOSr-"^ + Gr^ • r^^ . (12) 



Notice that the second term of S in ( 12 ) is the only new effect that the Schwarzschild geometry introduces 
on the previous result for the de Sitter case. The correction just comes from the Gaufi-Bonnet invariant 
part of i? in Q, which reduces to the square Riemann tensor G = RabcdR"'^^'^ = 12 r^/r^ for the 



Schwarzschild geometry, and induces the aforesaid r-dependence. From the expression of i? in (12 1, we 
see that in the limit of very large radius, r ^ r^-, the leading term of ([T]) that implements the relaxation 
mechanism is still the j3/B one, and this fixes /3 to be of order /3 ^ p\l^t ^ P\ -^o i ^s we have discussed 



in (10). Thus, in this large scale regime the Pq/R term plays no significant role. However, for local 
astrophysical scales the situation changes and then it is the 1/R term which takes over. Unfortunately, 
a detailed discussion of the SdS geometry is difficult. Even in the case where is set to zero, the field 
equations (|5| for the relaxation mechanism become quite complicated owing to the r-dependence of B. 
Albeit this case will be treated numerically in section |4] some qualitative considerations will be helpful 
to better grasp the behavior of the CC relaxation mechanism in the local astrophysical domain. In this 
simple context we can keep both the /3 and /J© terms. 

From the structure of the extra tensor E^b in ([3| , one can see that the R and B invariants appear in 
the denominator of Eab with some powers, schematically 

with ni,n2 > 0. Ignoring the details of /1.2 and ni^2 we may use the simple form 

°° 4 R^ B I6re^ ^ 108re-4 _^ g^2 . ^-6 ' ^^'^^ 

which should suffice to obtain some qualitative results. For r — > 00 we recover the de Sitter result in 
Eq. ([9]), which is an exact solution. For smaller values of r, the non-constant term 6r^ • starts to 
become more important and it will eventually dominate in (3/B, specifically for r < rd := (r^r^/lS)^/®. 
In this region the SdS metric is not a good approximate solution of ([S]) . 

Finally, for small r in the local astrophysical domain, we have 1/B cx r^ and the infiuence of this 
term wanes in comparison to the r-independent contribution 1/R ^ r^, which then takes over. Indeed, 
for small enough r, we find 



Eab^ [Pe^+f3:^j gab^ Pe^ 9ab for r « r,2 ( r^r^ ^ ) . (15) 



1/6 



In the local domain r <C rc2, the tensor Eab is approximately constant and the equations ([s]) are again 
satisfied for an appropriate dynamical choice of r^.. Thus, the SdS metric becomes an acceptable solution 



for dynamical relaxation in this region, too. From (151 it follows that the CC can be successfully relaxed 
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in the astrophysical domain provided is dynamically driven to fulfill the relation /3q r'^/8 + p\ = , 
which again in practice means p\ + (3q = 0, i.e. Eq. ([s]). In any case, we conclude 



1/30 1 ; 



(16) 



where = A4 . Compare (16l with the corresponding result in the cosmological domain, (lOl. 

Choosing Ai'^ = Ml Hq we get once more ~ H^^. In this case, however, for Mx 10^^~^^ GeV, the 
new (astrophysical) scale M.' must be located in the MeV range, which is a common scale for the SM of 
Particle Physics. It suggests a possible natural connection of the required order of magnitude value of the 
parameter Pq = —8p\r~^ ~ —p\Hq with conventional physics. Therefore, we could venture a possible 
natural explanation of why Ve is so large, also in the astrophysical domain: it might be only because 
Mx is much larger than the typical particle physics scales in the SM. If so, the hierarchy problem in 
Particle Physics could be linked to the hierarchy of vacuum energy scales in astrophysics and cosmology. 
However, this requires further studies. Far beyond the local astrophysical scales, we have the intermediate 
domain corresponding to the intergalactic and the intercluster distances, which is more complicated to 
analyze. See Sec. [4] for a preliminary study of that cosmological intermediate region. 



3 Studying the model in different gravity frames 

In this section, we investigate the existence and possible impact of extra d.o.f. in the relaxation model 
at small length scales, e.g. in the Solar system environment. In many theories of modified gravity new 
d.o.f. show up, which could mediate an extra force on matter coupled to gravity. A common way to 
identify the new d.o.f. from the modified gravity action is by considering its equivalent conformal scalar 
field representation. Using the arguments from Sec. [2] gravity in this regime can be described in good 
approximation by a modified F{R) theorv ^51 E51 126j . Indeed, the term f3/B involving the Gaufi- 
Bonnet scalar G is neglected and the the extra d.o.f. is just one scalar field. Furthermore, since our 
relaxation model is not just General Relativity with small corrections, it is convenient to discuss the 
situation in more detail. In particular, we will cross-check our results by considering a direct expansion 
of the F(R) theory around a constant background and then identifying the new d.o.f. 



3.1 Jordan frame 



The action ([ij is formulated in the Jordan frame, where matter is minimally coupled to gravity, i.e. to 
the metric gat- With (3 ~ we write it as a pure F{R) theory plus matter. 



S = 



-M 



2 / d^ 



-MlR~^ 
2 ^ R 



Pa 



Smat [9ab 



M'p / d''xV\g\FiR) + Smat[gabl 



(17) 
(18) 



where we have defined the reduced Planck mass Mp = l/V'8^ - IQi^GeV. As in Sec. § we use the 
SdS metric ( |11[ ) as a starting point, which allows fixing the parameter /3q. Then, the Einstein equations, 
emerging from 2SS/Sg'^^ — 0, are approximately given by 



— r 



Pa 



9ab 



0, 



(19) 



The equations (19 1 are approximate since 



where Tab follows from 2S{Smat)/Sg'''' and R = 12r- 
we keep the matter content and at the same time we use the de Sitter value of the curvature, which 
is justified inasmuch as we may use the de Sitter solution as an useful approximation to the current 
universe. Indeed, if we would strictly apply the condition (|6|, then the equations ( 19 1 boil down to (|5|, 
with Eab given by (15 1. We discuss the stability issue in Sec. |3.3[ Note that from Eq. (19 1 it follows 
that small changes in /Sq (or p\) yield only small changes in re and the solution does not change much. 
Thus, while in our case p)^ is canceled by a dynamical choice of Ve in which the matter sources do not 
participate in a significant way, in standard General Relativity the large value of p\ must be canceled ab 
initio by a very precisely chosen counterterm, otherwise the solution would change drastically, spoiling 
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completely the observed cosmological evolution. This observation displays the whole dilemma of the old 
CC problem, and points towards a further completion and improvement of the relaxation model. 

Next, following the standard procedure for converting an F{R) modified gravity model into a scalar- 
tensor theory, we introduce an auxiliary scalar field (/) by writing the action ( 18 1 as 

1 



AfJ, / d^x^\g\ [F{cj>) + F'{<t>){R - ^)] + ^matkb] 



S 



The original action can be recovered from the variational principle 6S/S(j) = 0, which yields 



F"{(j)) 7^ 0. For the model given in (18 1, we have 

/30 2pX 



^ M|,' 

where (3q — 2/3q/Mp. In the SdS solution, where R = <j> ^ 
reads 



(20) 



R if 



(21) 



^ 12rg ^, the diniensionless value of F'{(ps) 



Pa 



(22) 



with the critical energy density p° — 3Hq / {8tt Gn) — 3r~^Mp at present. Obviously, |F'((/)s)| ^ 1 is 
much larger in the relaxation scenario than in more common F(R) models [231 1241 126] . where General 
Relativity is slightly amended by a small correction. 



3.2 Einstein frame 

To obtain a standard Einstein-Hilbert term plus a scalar field responsible for the fifth force we apply the 
conformal transformation 

F'i^)gab = 9ab, (23) 

which takes us from the Jordan to the Einstein frame. 

Within the Einstein frame, endowed with the metric gab and corresponding curvature scalar R, the 
transformed action reads 



1 



S = -M 



,r2 



U I ^ -ab'^aF'WbF' 1 , 



yob 

F' 



The canonical kinetic term for the scalar ip results from the field redefinition — ^ (/s as follows: 

F'(0) = e 



exp 



2 ip 
3Mp 



Therefore, the action with canonically normalized fields is given by 



S = 



1 



-MAR 



1 



r''aa^a,^-^^(0F'-F(0)) 



9ab 



F' 



(24) 



(25) 



(26) 



from which the the scalar potential can be read off, and its value for the previous SdS solution be derived: 



Ml {ct>F'{cj,) - F(0)) 



Ml 



2{F') 



l\2 



{F'Y \ CP 



Ml 



F' 



,)' 



(27) 



with (j) — 4>s — 12rg ^ (associated to the value ips of ip in (25)), p\ ~ —^I'e ^ ^'^l F' (4>s), and moreover 
I3q = (12rJ^)^ _F'(0s). Obviously, V{ps) is the vacuum energy density in Einstein frame variables 
(coincident to the critical density in that frame, in the de Sitter approximation). 



3.3 Stability 



For studying the stability of the SdS solution we have to analyze the potential in (p7|), rewritten as 

Pa 1 „~-2x 



V{ip) ^MlJ/3Q\V~e 



1 



£2 £2 e 

4 



-2x 



(28) 
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was chosen such that U{0) = — |e^'/^ = —10. The circle is at the maximum and the square at the 
inflection point of U{x). To the right of the latter the scalar mass squared is positive. Note that both 
limits U{x ^ 0) and C/(a; — ?> oo) — > are valid only for the pure Pq/R model since the effects of the (i/B 
term will be important there. 



where from (22 1 we have |e^| ^ 1 around the SdS solution. Moreover, Pq — (12rg ^)^e^% in which 
_ F'{(j)s) according to (25l. Thus, within the same approximation, 



(29) 



Exactly at the SdS solution x — Xg we obtain V'{lPs) = and V"{lPs) — — 6r^^/e^= < 0, indicating that 
the SdS solution is actually a maximum, see Fig. [T] Noting that a scalar mass term in the Einstein frame 

{if'^/F') = a/I^I (F' m"^) (fi^ into the Jordan frame, it follows 



transforms as y'\g\ m — >■ (F'') 



that the mass squared of ip in the Jordan's frame is, 



m our case, m; 



LP,J 



The fact that it is 



negative is not a surprise, it was expected from the aforementioned fact that at the cosmological level the 
de Sitter solution is a maximum and hence unstable [145, its lifetime being of order Hq^ . However, when 
g2;/2 moves to larger values than ^e^=/^, then m? 



becomes positive and acquires physical meaning. From 
At the inflection point the vacuum energy 
y(e^/2 ^ 4g;r,/2) ^ (189/256) y((ps) is of the same order of magnitude as the de Sitter cosmological 



then on the scalar field (f behaves like decaying quintessence 

4 



constant. Thus the relaxation mechanism works also in the stable quintessence regime. 

It is interesting to reconfirm the above Jordan's frame calculation of the mass squared of the new 
gravitational d.o.f. by expanding around the constant background Ricci scalar Rq = 12r~^. Following 
the general procedure of [29] [30] j we obtain the vacuum Klein-Gordon equation in the Jordan frame. 



namely (□ 
squared: 



''Ri, 



Ri 



0, where the gravitational d.o.f. is denoted now by Ri{r), and m^^ its mass 
[F'(i?o) - RoF^iRo) - 3nF'(i?o)] 



3i^"(i?o) 



(30) 



and we used F{R) = -Pq/R, F'{R) = Pq/R^, F" {R) = -2/3© /i?^. We note that m^^ exactly coincides 
with TO^ J, which confirms our expectations. From this alternative approach, we most clearly see that 
the large value of Pq oc p\ does not enter the mass calculation, which means that the result does not 
depend neither on the size nor on the sign of p\. This is different from the previously considered jj!^ / R 
or R^lii^ models in the literature [ISKlllSSlliSllSTllMlIMlISO], where the mass is directly related to the 
parameter /i. 
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3.4 Coupling to matter 

Let us consider first the influence of matter on the evolution of tp, and leave for the next section the study 
of the role of ip in the field equation of matter. The corresponding Einstein's frame action reads 



1 



1 



1 



(31) 



where gab — gab/{F') is the (^-dependent Jordan frame metric, and £,nat is the matter Lagrangian density 
still in the original Jordan's frame variables and containing the density factor in the Jordan metric. 

First, we look at the Einstein equations following from 2SS/Sgab = 0, 



2 SCmatigab] 

Ad\ Sgab 



0, 



(32) 



where T^^ is the standard energy-momentum tensor of the canonically normalized scalar field including 
V{ip). The matter energy-momentum tensor reads 

2 SCniatigab] 2 (5£,„at [ffaf,] dgab 



2 SCraat[gab] 1 



rpab 

mat - yiii g^^^ gg^^ Q~^^ ^g^^ p, 

which depends on too. Next we discuss the equations of motion for p, given by 

— = -Uip -V {p) + — = 0, 

5(p dip 



(33) 



(34) 



where the variational derivative in the last term can be determined with the help of Eq. ( 33 ) 



5C^a.t[gab] _ 5Cn^at[gab] dgab 



5p 



Sgab dp 



p/rpab ~ ^'^ 
^ -^nuLtyab 



dp [F'Y 



V(iMp 



and in the last equality we used (25 1. Thus Eq. (34 1 can be cast as Dp + V'{p) 



rpab ~ 
^inatyab : 



b 7. 



(35) 



_ abl (V&Mp) = 0, 

where tp is seen to couple to the trace of The relation with the corresponding energy-momentum 



mat 



tensor in the Jordan frame follows from (33 1: 



{F'y 



r,"^ati ^'^d hence the 



_ ^ mat -VMM TS.it /F' 

traces in both frames are related by gabT^^at = (F') gabT^^ 
the physical (non-relativistic) matter density in the Jordan frame, we arrive at the effective field equation 
for p: 

np + V'{p) + ^ (36) 



Finally, since T^itgab = Pmat provides 



We may now compare the two source terms driving the evolution of p, viz. V'{p) and Pmat/ {MpF''^). If 
we look at the beginning of the quintessence regime {p > ps), where pmat ^ Pc 



r„ ' , the first source 

2\ 



is of order V'{p) ^ {Mpr^ '^)/F' ^ p^^/{Mp F'), whereas the second source is of order - p^J{MpF'^) and 
hence much smaller owing to the additional suppression factor of F' . This shows that the evolution of p 
is dominated by V{p). Hence, matter does not disturb the relaxation mechanism. This is not surprising 
because it is related to the dominance of the curvature by the relaxation mechanism. Indeed, if we take 



the trace of Eq. (32 1, we obtain MpR — gab + T^^^gab- From (27 1 we see that the first term on the 



righthand side ot' the trace is of the order V{ps) ~ (F') ^ durmg the slow-roll regime, whereas the 
second term is ^ (F')^^, which is once more smaller than V{ps) by a factor of {F')^^. 



3.5 Possible fifth forces with matter 

In the previous section we found that the evolution of the gravitational scalar p is not disturbed by 
matter. Conversely, now we wish to study if matter can be significantly affected by the presence of p. We 
will illustrate it by considering the effect on a matter scalar field ip with mass in the Jordan frame. 
The corresponding matter action reads 



t[gab,ip] 



^g^'dai^dti^-^ml^^ 



(37) 
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From the conformal transformation (23), and the associated relation g — {F')g , we switch into the 



Einstein frame. The complete action in this frame is given in Eq. (26 1, where in the present case the 
matter part reads 



mat — ^mat 



Sn 



gab 

F' 



1 1 



F 



1 1 



7 2r''9a^5,^- (^2^1^^ 



(38) 



The equation of motion for the matter scalar field ^ follows as usual from the variational principle 
SSinnt/Sip — 0. After some standard manipulations and partial integration, we arrive at 



5S^ 



1 



Si;. 



(39) 



The first term in the bracket at the integrand can be written as 



1 



— -+9 9i,ipda[ — ]=—a^P- —da'f, 



F' 



[F'Y dip 



and then using (251 we finally obtain the full equation of motion for the matter scalar 



(5bV) 



2 dgip 

3 Mp 



(40) 



(41) 



When the gravitational scalar tp is in a slow-roll regime we have Ida^p] ^ \/2 V{ips) ~ V^Mp ^ /^/F', 
where V{ps) was given in (27l. The maximal coupling with matter appears when da^p saturates the slow 



roll bound, in which case the equation of motion for ?/' in the Einstein frame can be cast in a simplified 
notation as follows: 



F'{p>) 



?A = 0. 



(42) 



Here, the mass m^^E = rn^/y/p^ corresponds to that of the matter field ip in the Einstein frame and 
contains the factor y/F' . We see from (42 1 that this factor is shared by the fifth force "friction" coupling 
of matter (i/i) with ip in that frame: T^^e ~ '''e^ I \/ F' {'f)- Since the dimensionless ratio T^^E/fn^^E = 
/m^ is independent of F' , the relative strength of the fifth force acting on the matter field ij} (measured 
by comparing it to its mass) is the same for all frames. Thus, we find F^^j = (F^.e/m^^^;) = 
for the effective fifth force strength in the Jordan frame. In other words, this coupling is essentially given 
by the present Hubble rate: Ki Hq. Even in the considered case, where the coupling of ip to matter 
is maximal, we obtained an extremely weak fifth force. Therefore, we must conclude that such coupling 
is completely unobservablej^ It is also interesting to mention the effect of ip on photon interactions with 
matter. Since the action for photons is invariant under conformal transformations of the metric, they are 
not subject to explicit fifth forces mediated by ip. However, they feel the metric deviation in our modified 
gravity model with respect to standard General Relativity. This effect might be detected by gravitational 
lensing experiments depending on the strength of the deviation. We briefly address this possibility in the 
next section. 



4 Numerical results for the large scale astrophysical domain 

The large scale cosmological domain characterized by the FLRW metric was studied in detail in |14| [T5]. 
and we saw it is dominated by the (3/3 term in ([ij. We have argued in Sec. [2] that the (3/B term in the 
action is also the dominant term at large length scales for the metric ([?]), which is more in accordance 
with the astrophysical environment. Unfortunately, due to the occurrence of the Gaufi-Bonnet invariant 
in ^ it is difficult to find exact solutions on a spherically symmetric background. This is quite obvious 
from the complicated structure of the field equations ([2|-([3]), especially when the curvature invariants are 
non-constant. We will show here some numerical solutions of these equations in the vacuum case and we 
briefly discuss the deviation from the standard SdS solution. The latter metric, as given by equations 
([7| and (111, will be used for setting the initial conditions for large radius r, i.e. the situation for which 
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Figure 2: Top: numerical solutions for the functions a{r) and b{r) in the spherically symmetric metric ([T]) 
for the pure P/B model without matter. The initial conditions have been fixed close to the higher zero 
of the corresponding SdS solution (111, which is plotted, too. Bottom: relative differences 
the numerical functions a{r) and b(r) 
of both functions. 



-b(r) 



of 



a(r)+6(r) 

Note that the numerical error dominates around the higher zeros 



the /3-terms dominate over the i3q ones. This procedure ensures a smooth transition to the cosmological 
de Sitter result. 

Specifically, to obtain the pure effects of the /3/B term, we set (3q to zero in ([3|, and the initial 



conditions for the metric are set very close but below the higher zero of the function a{r) in Eq. (111. For 
small values of the Schwarzschild radius parameter <C this is approximately aX r — r^. As we know 
already, r^. = corresponds to the de Sitter solution, which is exact also for the ji/B model. Therefore, we 



■'One can show that if we would also transform the matter field tp into the Einstein frame, i/i — ^ \rF' -ij), then the tiny 
fifth force could be fully absorbed into a very small, unobservable, correction to the matter field mass. 
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consider only solutions with non-zero r^. Some examples are shown in Fig.|2] where the two functions a(r) 
and 6(r) in the metric ([t]) are displayed together with the standard SdS result. Let us note that in the case 
of small perturbations, one usually writes a{r) = l + 2'i/{r) and b{r) = 1 + 2 $(r) where ^{r) and $(r) are 
the two gravitational potentials associated to the metric. While in the standard SdS case they are equal, 
in the present case we look for deviations a{r) — b{r) = 2[\['(r) — $(r)] that could provide information of 
the modified gravity model. As expected, these deviations grow with increasing r^. While b{r) follows 
qualitatively the standard result, the function a{r) stays at larger values especially for small r, which 
we identify as an effect of the r-dependence in B. The relative difference between both metric functions 
is plotted also in Fig. [2] (bottom), it is of order one around the lower zero of b{r). However, we expect 
that this singular point must lie within the matter distribution, whose metric has to be matched to the 
vacuum solution in Fig. [2] at a larger radius r. Whether this matching is possible and how the inner 
solution for the metric behaves, depends strongly on the properties of the matter distribution. We keep 
this question open for a future investigation. At the moment we cannot exclude the possibility that at 
intermediate (intergalactic) scales some non-trivial corrections to standard expectations might creep in. 
We provide some discussion of this possibility in the next section. 



5 Possible gravity modifications at large distances 

Let us return for a moment to the discussions made in Sect. [21 We have seen that while both at local 
astrophysical scales and at very large cosmological domains the relaxation mechanism works perfectly 
well and can be studied analytically, at intermediate scales our knowledge is more limited and we need 
to resort to the numerical analysis within some approximations, as we have seen in the previous section. 
The characteristic distance entering the potential deviations is 

r>lf) , (43) 

where we know that /3 and /3q are fixed in order of magnitude and (same) sign by the condition that 
the relaxation mechanism works in the cosmological and astrophysical domains, respectively. These 



conditions entail the relations \/3\ ^ Pa^o ^^^^ l/^ol ^ Pa^o- Hence it follows from (43 1 that 

(44) 



1/6 xi/3 /G^M,^^/^ 



where Ms is the mass of the galaxy or cluster of galaxies sourcing the gravitational field. We may 
compare this length scale with the one at which one expects deviations to the Newtonian gravity at large 
distances in ordinary extended gravity theories, i.e. with a late-time modification of the gravitational 
interaction These effects have been specifically addressed in Ref. by 

considering a Lagrangian modification of the form SC = Mp /i'^'^+^/Q" with Q = R°'^'^'^ Rabcd the square 
Riemann tensor. For n = 1 this is not far away from the behavior of the gravity modification in our 
relaxation model ([T]) at large distances, except that in our case we have both a gravity modification and 
the large CC in mutual interplay. It is however illustrative to compare the two kind of models for n = 
for which the the characteristic scale of the induced modification of Newtonian gravity is 



9\ 1/5 

{GNMsf \ 
HI i 



(45) 



after setting ^ Hq, which is the extremely small value that is necessary to choose for the /i-parameter 
in order to explain the late-time cosmic acceleration in this kind of models [31 . The ratio between the 
two characteristic length scales is 



(GatM.Ho)'^'' =0(0.1). (46) 



^The CC relaxation model has also been analyzed at the cosmological level in more general cases in which we have an 
arbitrary power of the denominat or and also when a power of the Ricci scalar R is included in the numerator of the 
gravity modification in ([TJ, see |14II15I for details. However, here we just analyzed the canonical, i.e. the simplest realization 
of these models at the astrophysical level, and therefore we should naturally compare also with the simplest case n = 1 of 
Ref. |3T|. 
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For the numerical evaluation of this ratio we have used the source mass of a relatively large galaxy (say, 
Ms ~ 10^^ Mq, i.e. with a number of solar masses similar to our galaxy), or the mass of a large cluster of 
similar galaxies. In both cases the two length scales do not differ in more than one order of magnitude, 
so they are indeed very close. In other words, in both classes of models (relaxation model and ordinary 
late-time extended gravity) we expect deviations from Newtonian gravity at similar large scales in the 
100 kpc-Mpc range, as it follows from the above formulae. Let us remark, however, that in the relaxation 
model we have a double bonus which is absent in the ordinary case, to wit: i) we do not need to use 
extremely small scales as /i ~ Hq ~ 10"'^^ eV but rather scales M. ~ meV and M.' ~ MeV both lying in 
the natural SM range (cf. Sect. [2]); and ii) the local astrophysical domain is protected from unwanted 
vacuum effects even in the presence of the large p\. Indeed, the typical value of Tc emerging from (44 1 
is roughly from 10 to 100 kpc in the case of a galaxy, whereas for a large cluster it may reach the few 
Mpc level. As these length scales are at the same time the characteristic sizes of a galaxy and a cluster 
of galaxies, respectively, it follows that the local astrophysical domains remain well protected from large 
vacuum effects in our relaxation models. Finally, in the intermediate region beyond these distances (say 
up to '-^ 10 — 100 Mpc) we cannot study the problem in a simple way, as the numerical analysis of the 
previous section has shown, because the SdS metric is no longer a good approximate solution. Only after 
attaining very large distances (of order of a few hundred Mpc at least) we retrieve again the relaxation 
mechanism in the cosmological domain [2] with all its phenomenological success [T5] . 



6 Conclusions 

In this letter, we have performed a first investigation of the modified gravity implementation of the CC 
relaxation mechanism in a static and spherically symmetric background, which serves as an approximation 
for the study of the mechanism at length scales characteristic of the astrophysical domain. We found the 
mechanism to be working in the sense that the space-time curvature is not dominated by the initial large 
vacuum energy density p\. Instead, a small effective dark energy density 0{Mpr~^) ^ Mp ~ is 
found, similar to our earlier results in a cosmological background |14| . Moreover, we studied the influence 
of fifth forces on massive objects, which turned out to be long ranged, but extremely weak and therefore 
inconspicuous. This shows that the idea of CC relaxation may be useful in more general contexts than 
just in cosmology. Given the unconventional starting point, one could not naively expect this local 
astrophysical behavior from the beginning. Interestingly, we also found that at very large scales of order 
of 10~^ Hq^ the metric may deviate significantly from the standard SdS solution, thus opening up a 
possible way to detect this mechanism via gravitational lensing and other effects. 

Possible observable implications at these scales have already been proposed in other modified gravity 
theories [31]. In the light of these studies one can foresee deviations from the standard Newtonian/GR 
behavior at length scales of 100 kpc-Mpc. We cannot exclude, for instance, the possible implications of 
these deviations e.g. by changing the required abundances of dark matter in the intergalactic domains in 
order to explain the DM problem in the large. However, the investigation of these potentially relevant 
issues goes far beyond the content of the present Letter. Their more precise description is therefore an 
important future challenge in the CC relaxation approach. 

Finally, some theoretical issues remain open for further study; in particular, the fact that the relaxation 
mechanism exerts such a dominating influence. In this sense, further insight should go into the direction 
of finding ways to complete this model and understanding the role of the vacuum energy versus the matter 
sources. Quite in contrast to the present approach, the usual modifications of gravity disclaim ab initio any 
understanding of how to cope with the huge vacuum energy injected in the universe by QFT/string theory, 
and just focus on late-time simulations of the measured effect. This way is therefore no more satisfactory. 
A link between the two points of view is still missing, and hence more work is required to understand 
the ultimate interplay between gravitation, vacuum energy and matter in our cosmos. We hope that our 
approach may offer a new perspective for an eventual solution of this difficult problem. 
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